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Abstract 

The linear response of a nonrelativistic superfluid baryon system on an external weak field is 
investigated while taking into account of the Fermi-liquid interactions. We generalize the theory 
developed by Leggett for a superfluid Fermi-liquid at finite temperature to the case of timelike 
momentum transfer typical of the problem of neutrino emission from neutron stars. A space-like 
kinematics is also analysed for completeness and compared with known results. 

We use the obtained response functions to derive the neutrino energy losses caused by recom- 
bination of broken pairs in the electrically neutral superfluid baryon matter. We find that the 
dominant neutrino radiation occurs through the axial-vector neutral currents. The emissivity is 
found to be of the same order as in the BCS approximation, but the details of its temperature 
dependence are modified by the Fermi-liquid interactions. 

The role of electromagnetic correlations in the pairing case of protons interacting with the elec- 
tron background is discussed in the conclusion. 
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I. INTRODUCTION 



Thermal excitations in superfluid baryon matter of neutron stars, in the form of broken 
Cooper pairs, can recombine into the condensate by emitting neutrino pairs via neutral weak 
currents. This process was suggested [l| many years ago as an efficient mechanism for cooling 
of neutron stars in some ranges of ternperature and/or matter density. The interest in this 
process has been recently revived in connection with the fact that the existing theory 

of thermal neutrino radiation from superfluid neutron matter leads to a rapid cooling of the 
neutron star crust, which is in dramatic discrepancy with the observed data of superbursts 
Q], 0. It was realized that a better understanding of an efficiency of the neutrino emission 
in the pair recombination is necessary to explain modern observations. 

The relevant input for calculation of neutrino energy losses from the medium is the 
imaginary part of the retarded weak-polarization tensor intimately connected with the au- 
tocorrelation function of weak currents in the medium. Though the theoretical investigation 
of the autocorrelation functions of strong-interacting superfluid fermions was started more 
than four decades ago the complete theory of the problem does not yet exist. Leggett's the- 
ory of a superfluid Fermi liquid is limited to the case when both the transferred energy 
and momentum are small compared to the superfluid energy gap, i.e., cj, q <^ A. This theory 
cannot be applied to calculations of neutrino energy losses because, in this case, we need the 
medium response onto an external neutrino fleld in the time-like kinematic domain, u > q, 
and u > 2 A, as required by the total energy u = ui + UJ2 and momentum q = q^^ + q2 of 
escaping neutrino pair. 



The well-known Larkin-Migdal theory 



temperature. 
5|, p, where 



,10] is restricted to the case of zero 
Recently, the calculation of the neutrino energy losses was undertaken in Refs. 
the imaginary part of the autocorrelation functions was calculated for a superfluid neutron 
matter at zero temperature. This approach is apparently inconsistent because the imaginary 
part of retarded polarization functions substantially depends on the temperature [see Eqs. 
( l82l) . ( l88l) . (11121) . and (11131) of this work]. One more inconsistency of the work [q!] is, that 
the temporal component of the axial-vector current cannot be discarded, as it is done by the 
authors. This relativistic correction contributes to the neutrino energy losses of the same 
order as the spin-density fluctuations, i.e. oc Vp. This was pointed out for the flrst time in 
Ref. III]. Below, we will return to the discussion of these works and compare our result 
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with that obtained in Refs. [sl, and in some earher works. 

The appropriate, temperature-dependent approach is developed in Ref. {4], where the 
mean-field BCS approximation is used to calculate the superfiuid response in the vector 
channel. To include the Fermi-liquid effects discarded in the BCS approximation, in this 
paper, we first generalize Leggett's theory to the case of arbitrary momentum transfer. We 
evaluate the weak-interaction effective vertices and the autocorrelation functions while tak- 
ing into account strong residual particle-hole interactions. To obtain a solution of Leggett's 
equations in reasonably simple form, we approximate the particle-hole interactions by its 
first two harmonics with the aid of the usual Landau parameters. Within these constraints 
we obtain the general expression for the autocorrelation functions and then focus on the 
superfiuid response in the time-like kinematic domain. We investigate both the vector chan- 
nel and the axial channel of weak interactions to evaluate the rate of neutrino energy loss 
through neutral weak currents caused by recombination of electrically neutral baryons. 

The role of electromagnetic correlations in the pairing case of charged baryons interacting 
with the electron background deserves a separate consideration. The quantum transitions 
of charged quasiparticles can excite background electrons, thus inducing the neutrino-pair 



emission by the electron plasma 



12j, 131]. In summary, we briefiy discuss this problem in 



the light of modern theory to understand whether the plasma effects can lead to noticeable 
neutrino energy losses through the vector channel. 

The paper is organized as follows. The next section contains some preliminary notes 
and outlines some important properties of Green's functions and the one-loop integrals 
used below. In Sec. Ill, we discuss the set of equations derived by Leggett for calculation 
of correlation functions of a superfiuid Fermi liquid at finite temperature. In Sec. IV, 
we consider the superfiuid response in the vector channel. Because of the conservation 
of the vector weak current it is sufficient to consider only the longitudinal and transverse 
autocorrelation functions. The correlation functions in the axial channel are evaluated in Sec. 
V. As an application of our findings, in Sec. VI, we evaluate neutrino energy losses through 
neutral weak currents caused by the pair recombination in superfiuid neutron matter. Some 
numerical estimates of the neutrino energy losses are represented in Sec. VII. Section VIII 
contains a short summary of our findings and the conclusion. 

In this work we use the standard model of weak interactions, the system of units h = c = 1, 
and the Boltzmann constant ks = 1- 
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II. PRELIMINARY NOTES AND NOTATION 



In our analysis, we will use the fact that the Fermi-liquid interactions do not interfere 
with the pairing phenomenon if approximate hole-particle symmetry is maintained in the 
system; i.e. the Fermi-liquid interactions remain unchanged upon pairing. According to 
Landau's theory, near the Fermi surface, p ~ p'~ p^, the Fermi-liquid interactions can be 
reduced to the interactions in the particle-hole channel. We will assume that the effective 
interaction amplitude is the function of the angle between incoming momenta p and p' and 
can be parametrized as the sum of the scalar and exchange terms 

aVf- (nn') = / (nn') + g (nn') Y.M[- (1) 

Here and below, p = pirM*/7r^ is the density of states near the Fermi surface; n = p/p and 
n' = p'/ p' are the unit vectors specifying directions of incoming momenta, a ~ 1 is a usual 
Green's function renormalization constant independent of u, q, and T, and o"j {i = 1, 2, 3) 
stand for Pauli spin matrices. The pairing interaction, irreducible in the channel of two 
quasiparticles, is renormalized in the same manner 

a'pt^ (nn') = (nn') + F^ (nn') ZMl- (2) 

We will consider the case when the pairing occurs only between two quasiparticles with the 
total spin S = 0. Then the irreducible pairing amplitude is to be taken as the singlet, 

a^pf ^ (nn') ^ F^ (nn') = F^ (nn') - 3F^ (nn') . (3) 

Since the baryonic component of stellar matter is in thermal equilibrium at some tem- 
perature T, we adopt the Matsubara Green's functions for the description of the superfluid 
condensate and for evaluation of the polarization tensor. In the case of ^5*0 pairing, near the 



Fermi surface, these are given by: 



G (Pn, P) = a J\ l" , Gh (Pn, P) = a 'fr, ^ (P- P) = ^ ^2^ 

where p„ = tt (2n + 1) T with n = 0, ±1, ±2, ... is the fermionic Matsubara frequency. In the 
above equation, G and Gh represent the propagators of a particle and of a hole, respectively, 
and F is the anomalous propagator, i.e. the amplitude of the quasiparticle transition into a 
hole and a correlated pair. For the inverse process: (p„, p) = F {pn, p). 



We use the momentum representation and the following notation 

- _ 'Pf _ \ (r\ 

where M* = pp/Vp is the effective mass of a quasiparticle, and the energy of a quasiparticle 
is 



Er, = ^el + AHT). (6) 

We designate as Lx,x q; p) the analytical continuation onto the upper-half plane of 
complex variable u of the following Matsubara sums: 

LxX' {uJm, P+|; P-|) =T^X (^Pn + UJm, P+|) ^' (Pn, ) ' C^) 

Pn 

where X,X' G G,F,G^. 

In Leggett's equations, which we are going to exploit, the spin dependence is already 
taken into account, and Xlpo- everywhere replaced by 2^p. It is convenient to divide 
the integration over the momentum space into the integration over the solid angle and over 
the energy according 
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Pf'^l <fep-- (S) 



(27r)' J 47r 

and operate with integrals over the quasiparticle energy: 

Ixx' iuj,qcos9,T) = J dSpLxx' (uj,p+^,p-^^ . (9) 

These are functions of u, q, and cos^^ = nnq, which is the polar angle of the direction of the 
momentum p = pn relative to the direction of nq = q/q as the z axis. 

The functions 2xx' possess the following properties, which can be derived by a straight- 
forward calculation |9|: 

^GF = —^FG, ^FGh = ~^GhF, (10) 
^GhF + IfG = -^FF, (11) 
'^GhF — IfG = —^^FF, (12) 

- {Igg, + Tff) =Ao+ ^-^^^^^Tff. (13) 
Here v =Vpn, and the quantity Aq = — (IgGh + ^FF)q=ow=o satisfies the gap equation 

1 - T^Ao = 0, (14) 



where Fq is the zeroth harmonic of the singlet pairing amphtude ([3]). 

The key role in the medium response theory belongs to the functions defined 
following combinations of the above loop integrals: 

A (cu, qn) = a~'^Xpp, 

K {uj, qn) = a"^ Q (Tgg + ^g'-g'-) + ^ff 

X qn) = a'^^ {Igg - Ig^g^) ■ 
These can be derived in the following form: 

A2 r de 



^ = I /I ^ " + + ^-^ ■ 

To shorten the expressions, we use the following notation: 
and 

1 f , E+ . E\ 

$4- = ^ ^ tanh — - ± tanh — - . 

{cu + iOf -{E+±E_f \ 2T 2TJ 

It is straightforward to verify that 

A {u, qn) = A {ui, — qn) , k, {uj, qn) = k {uj, — qn) , 
and that the functions n {uj, qn) and x Qn) are not independent because 

UJK = qvx- 

III. LEGGETT'S FINITE-TEMPERATURE FORMALISM 

The two-particle autocorrelation function is defined as 



{Uj,q) = 2^ e (P,C^) ( («p+q/2,.«p-q/2,<x ■ V-q/2,a' V+q/2,<x' ) ) ^ (p',C^') , 
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where ^ (p,cr) is a three-point vertex responsible for the interaction of a free particle with the 
weak external field. It is some function of the momentum p and spin variables a; ((^4 : B))^ 
is the Fourier transform of a retarded two-particle Green's function. 

The analytic form of the autocorrelation function can be immediately written, if we know 
the effective (full) three-point vertices defined via the linear correction to the quasiparticle 
self-energy S^^^ (V) in the external field V (see, e.g., Ref. {l5|): 

T = e(P,o^) + ^. (26) 

Near the Fermi surface, these vertices can be treated as functions of transferred energy and 
momentum, q = {u, q), and the direction of nucleon motion n. 

In superfiuids, we have to distinguish the vertices of a particle and a hole, which are 
related as C,h (p,o") = ^ (— P, —o'). Since there are two possible cases, ^ (— p, —cr) = ±^ (p,c"), 
it is convenient to consider the "even" and "odd" bare vertices 

e±(n) = i(e(p,^)±e(-p,-^)). (27) 

We denote as 

T±(n) = i(T(p,a)±T(-p,-a)) (28) 

the corresponding full vertices taking into account the polarization of superfluid Fermi liquid 
under the influence of the external field. 

In Eq. f l26l) . the quasiparticle self-energy consists of the normal part and the anomalous 
part caused by the pair condensation. In the case of ^Sq pairing, the anomalous self-energy 
is sensitive only to the longitudinal vector fields, because the only kind of motion possible for 
the condensate is potential flow, i.e., a density fluctuation [16]. Therefore for the longitudinal 
currents, along with the ordinary vertices 7±, it is necessary to consider the anomalous vertex 
T, responsible for excitations of the condensate. 

As was derived by Leggett (see Eqs. (22) and (23) of Ref. [qJ), the longitudinal effec- 
tive vertices T±,T are to be found from the following equations (we omit for brevity the 
dependence of functions on u and q) : 

f(n)- I —r^(nn')Aof(n')- I ^F^nn') ^^^^^fc^A fn') T fn') 



47r ' ' ^ ' ' 7 47r ' ' 2A2 

. ^ (nnO (nO T (n') - ^ / — 

Att ^ ' A ^ ' ^ ' A An 



+ I ^F- (nn') ^A (n') T (n') [ ^F- (nn') A (n') (n') = 0, (29) 
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T_(n) + I ^/(nn') (n') t (n') - J ^/ (nn') (n') T_ (n') 

+ [^f (nn') -^K (n') (n') = ^ (n) , (30) 
J 4:71 qv' 

r+ W - I / ^/ (nn') A (n') f (n') + / ^/ (nn') (n') r_ (n') 

- I ^/ (nn') {k (n') - 2A (n')) (n') = (n) . (31) 
In Eq. fl29|) . the irreducible pairing amplitude is to be taken as the singlet, as given by Eq. 



Once the effective vertices are calculated, the two-particle autocorrelation function can 
be immediately found using the expressions: 



, dn 



^A (n) t (n) + —K (n) T (n) + (n (n) - 2A (n)) T+ (n) 
A qv 



(32) 



if ^ = and 



dn 



— \ (n) t (n) + ft (n) T_ (n) - —ft (n) T+ (n) 
A qv 



(33) 



One can easily verify that these equations represent a generalization for the case of finite 
temperatures of the Larkin-Migdal [10| equations derived in the ladder approximation for 
the vertices modified by strong interactions in a superfluid Fermi liquid. 

Unless we are dealing with a spin-independent longitudinal field only fluctuations of 
the normal component contribute to the polarization. The corresponding effective vertices 
should be found from the equations jol: 



71 (n) + 



r (n) - 
dn' 



dn' 

Atx 



fnn'l 



ftfn')r_fn') 



qv' 



ft fn') fn') 



(n) 



(34) 



47r 



fnn') 



UJ 



qv 



;ft(n')r_(n')-(ft(n')-2A(n'))r+(n') 



^+(n), (35) 



which represent Dyson's equations ideally summing the particle-hole irreducible diagrams 
in the ladder approximation. In these equations, the spin dependence is already taken into 
account, so ^ is to be taken as a function of only p, i.e. ^ (p) = ^ (p,o") = ^ (p, —a) for 
S" = 0, and ^ (p) = ^ (p,cr) = — ^ (p, — cr) for S* = 1. The c number refers to the 
usual Landau " quasiparticle-irreducible" scattering amplitude F"^ (n, n') as defined in the 
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normal phase; it is to be taken as the spin-independent or spin-dependent part according to 
e(p,<j)=±e(p, -a). 

In this way one may calculate the spin, transverse-current, and helicity-current autocor- 
relation functions, which are given by the expressions: 



/dn , , 



qv 



K (n) r_ (n) + {k (n) - 2A (n)) r+ (n) 



\ii = C+, and 



K (n) r_ (n) K (n) 71 (n) 

qv 



(36) 



(37) 



We arc now in a position to evaluate the autocorrelation functions necessary for calculat- 
ing the energy losses from a hot superfiuid baryon matter. We consider the medium response 
in the vector and axial-vector channels which are responsible for the neutrino interactions 
with the medium through neutral weak currents. 



IV. VECTOR CHANNEL 



Vector current of a quasiparticle J'* is a vector in Dirac space (/i = 0, 1, 2, 3). The corre- 
sponding polarization tensor Ily^ (a;, q) must obey the current conservation conditions: 

n^^'^ (a;, q) = 0, q^Xl^ (a;, q) = 0. (38) 

These equations imply that the polarization tensor can be represented as the sum of longi- 
tudinal (with respect to q) and transverse components 

(^, q) = (^, q) {l, ^n^^ ' {l, ^n^^ + (^, q) g^' {b^^ - n\n^^ (39) 

In this expansion, the longitudinal and transverse polarization functions are defined as 

ni = n°°, nT = ^((5'^-nX)tf^\ (4o) 

The transverse polarization function can be conveniently evaluated in the reference frame 
where the z axis is pointed along the transferred momentum, so that nq = (0, 0, 1). Then 

Ht (c^, q) = ^ (n^'i (u;, q) + tf'^ (c,, q)) . (41) 

Thus we actually need to calculate only the temporal and transverse components of the 
effective vertices. 



A. Longitudinal polarization 



The vector current of a free particle has the nonrelativistic form 

= (1, V) , (42) 
where v = p/M is the particle velocity. In this case we find 

e=c = 1, e = 0, (43) 

e+ = 0, e = t = v\ (44) 

Then the longitudinal polarization, II^, = Ki {uj, q) ,can be calculated with the aid of Eqs. 



(I29D-(|32D with ^+ = 1 and ^_ = 0. 

Before proceeding to the detailed solution of these equations, let us note that apart 
from the ground state, Eq. fl29l) allows for excitations of the bound pairs with the orbital 
momentum / > 0, if these exist. We will consider the simplest case of ^Sq pairing, assuming 
that the only possible bound state of the pair corresponds to the zero angular momentum I. 
This allows us to consider only the zeroth harmonic of the pairing interaction. In this case 
the anomalous vertex is independent of the quasiparticle momentum and the use of the gap 
equation f|T^ allows us to recast Eq. fl29|) as follows: 

Using Eq. (H3!l . we obtain Eqs. (!30|) and (!3T]) in the form 

T_ (n) + l^f (nn') (n') t (n') " / ^/ (nn') (n') T_ (n') 

+ f^f (nn') -^K (n') (n') = 0, (46) 

(n) - ^ / ^/ (nn') A (n') t (n') + / ^/ (nn') (n') T (n') 

- 1 ^/ (nn') {k (n') - 2A (n')) (n') = 1. (47) 

The vertex equations can be further simplified in various assumptions about the amplitude 
of the particle-hole interaction ([T]), which can be expanded in the Legendre polynomials, 
according to 

/ (nn') = J2 fiPi (nn') • (48) 
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We consider a simplified model with = for / > 2, when the interaction function is given 

as 

/(nn') = /o + /inn'. (49) 
Solution to the set of Eqs. fl45l) - fH7l) can be written with the aid of the following notation: 

a(^,q,T) = J ^X{n), 7(u;,q,T) = J ^A(n)cos^^, 

7] {u, q, T) = J (n) , /3 (u, q,T) = J (n) cos^ 9, (50) 



where 



s^a — 7 s^a — 7 



After some algebra, we find: 



(51) 



^+ l-f,(l + f,(s^Q-P))Q-f,P^ (^2) 
^ ^ sfiQcosO 

1 -/o(i + A (s2g-p))g-/iP' ^ ^ 

^^2- gMc:^(l-/i/?)+7r?/i) ( . 

uj{s'a-^)[l-h{l + Ms'Q-P))Q-hPy ^ ^ 

A short calculation of the right-hand side of Eq. (!32l) with ^+ = 1 gives the simple result 

"^'"■"■^' = ^-/o(i + /,(.^'g-P))Q-/.P '''' 

L 5C5 limit. 

Notice that the autocorrelation function of the density fluctuations has already been 
calculated in various limits. Let us take, for example, the BCS limit by setting /o = /i = 0. 
We then obtain 



This expression is in agreement with Eq. (37) of Ref. {4] if we take into account the relations 
77 — 2a = Aqo, t^tt = — AAq, and q^V^7 = AgjAj connecting our notations and those of Ref. 
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2. Limit uj, qVp <C A, T > 0. 



In this limiting case, from Eqs. (fTSl) . (fT9|) . and (150!) we find (see also Ref. 

1 1 /• rfn /"^ , (cos2 ^ - s^) e^lE'^ dn 



2 2 J Att Jo -{cos'^ 6) eyE^dE' 

1 1 f dn ^ cos^ ^ (cos2 ^ - s^) ^VE^ dn 



where 



Then Eq. (l55|l gives 



6 2 7 47r 7o s2 - (cos2 0) £2/^2 dE' 

rfn /-"^ , (cos2 0)e2/^2 
^ / de- 



Air Jq s2 - {cos^ 6) e^E^dE'' 



/? ~ 27 + s^r] - ^, 



(in 1 , -2 E 
-TT^ = cosh — -. 
dE 2T 2T 



1 - [/o + /is2/ (1 + /i/3)] Q (s) 
in agreement with the result of Leggett . 



3. Limit to, qVp <. A, T = 0. 



In the case T = 0, Eqs. (|57D-(|60D give 



1 1 

a = 2' 7 = g,^ = = 0, 



so 

We then obtain 



u ( T/ A T m p(l + /i/3)q%V3. 

IIl (cu, qVp < A, T = 0) - 



a;2-(l + /o)(l + A/3)q2Vj/3- 



in agreement with the results obtained in Ref. lOj. 
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4- Time-like momentum transfer, < T < Tc 

We are interested in the case of time- like momentum transfer, q< u, and uj > 2 A taking 
place in kinematics of the neutrino-pair emission. Then we deal with the case qVp <^ u, 
i.e., u = s^^ <^ 1. In this limit, we have 

^-67 = -Rea, j3 ^ rj u^a, (67) 

Using this fact, we find the functions Q and P in the forms 

Q = r] + 2u^-i, (68) 
s^Q- p = s'^r^- [3 + 27 (69) 

The real and imaginary parts of the functions can be obtained from Eqs. (fTSl) . ([T9l) . and 



(1501) . The real part can be evaluated to the lowest accuracy. We find: 

/oo ^2 

1 r de E 

Re„=!f!fl + 2J>r^^^ta„h^V (72) 

' 3 V J_^Euj^-AE^ 2T J ' ^ ^ 

Re^4(l + 2p£|-^ta„h|), (73) 

where the symbol V means principal value of the integral. In deriving the last two equalities 
we used the identity 

- — -qrrtanh — - / d6—— = l. (74) 
2j_^EE^ 2T 2j,^ E^dE ^ > 

Within a time-like momentum transfer and uj > 2 A, the imaginary part of the function 
arises because of the pole at = -Ep+q + Ep. We calculate the imaginary contributions up 
to the higher accuracy and find 

Im a = TT — P) (u — 2A) tanh — - 

uVoo^ - 4A2 ^ ' AT 

1 2 fuj^ + AA^ cu2-4A2 u \ 



M 1 + - -^h- - J ) , (75) 
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Im7 = ^ Pi (to - 2A) tanh — 

3 Vc^2 _ 4^2 ^ ' 4T 

Imr? = =B (tu — 2A) tanh —— 

6 2 /cu2 + 2A2 Cj2_4^2 ^ _^ ^ 



X 1 + — ^ ^ cosh"' — , (77) 

^ 5 V - 4A2 32T2 4T ' ' ' ^ ^ 



Im/? = -— =9 (a; - 2A) tanh — 

5 ujy/u^ - 4A2 ^ ^ 4T 



25 2 ^cu2 + 4A2 cu2-4A2 , _2 \ 



M 1 + [ijrrj^. - ^5t^ -^h- - J ) , (78) 



where (x) is the ordinary Heaviside step function. 
We also find: 



^ 2n u^A^Q {u - 2A) , a; 

g = % , ^ ^ tanh— , (79) 

^35 ujVuj^ - 4A2 4r' ^ ^ 

2^ „ 1 .Stt uWe(cj-2A) , u; 

3 14 (co'2 — 4A2) vco'2 — 4A2 4T 

^2 ^ 4^2 ^2 _ 4^2 _2 ^ 



and 



P = — [l + -V / ^ ^ tanh — 

-I , ^ ' tanh — ■ (81) 

5 uj^/uj^ - 4A2 4T ^ ^ 

Having these formulas at hand, we can evaluate the real and imaginary parts of the 

longitudinal polarization function fl55|) . After a little algebra, we obtain 

n,.(u;,q,T) = plvj(^l + i/i^^ 

Stt.mA 1 ,\^qWe(cJ-2A) , uj 

As one can see from this expression the spherical harmonic of the pairing interaction does 
not affect the longitudinal polarization in the high-frequency limit uj ^ qVp. If we set 
0, this expression reproduces the result of the BCS approximation [see Eq. (48) in Ref. 



m 
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B. Transverse polarization 



As explained above, the transverse field does not affect the anomalous self-energy of a 
quasiparticle. Therefore the transverse-current autocorrelation function 

Kt {uj, q) = ^ {Ki-=vi (w, q) + K^-=V2 q)) • (83) 

can be evaluated with the aid of Eqs. flMl) . (I35l) . and fl371) with ,^_|_ = and = v]_, where 
vj_ = (f sin 6' cos <y9, f sin ^ sin 0). The particle-hole interaction (H9l) can be written as 

/o + /inn' = fo + fi (cos 6 cos 6' + sin ^ sin 6' cos (v? - <^')) • (§4) 

The sets of equations for different z = (1, 2) are decoupled, and we find: 



and 



In the case q < cj, and c<j > 2 A, using Eqs. (172!) . (!73|) . (1771) . and (ITHj) . we find 
7]- p = —u^ {1 + 2V — tanh — 



15 V .Loo Euj^- 4^2 2T 



-2— — tanh—. (87) 

15 ujVu^ - 4A2 4T ^ ^ 

Up to accuracy Vp from Eq. ( l86l) . we obtain 

(u;, q) = —pV'^ (l + 2V r - tanh — ^ 

pV^^ tanh — . 88) 

». „.... : .St:,,.!., a ». .. .... 

frequency limit, u <^ 1, the first two harmonics of the particle- hole interaction do not affect 
the transverse polarization of the medium. 



V. AXIAL CHANNEL 

Since only the normal component contributes to the spin fluctuations, the axial effective 
vertices should be found from Eqs. and (135|) . and the corresponding correlation functions 
are given by Eqs. (!36|) and (l37|l . We now focus on this calculation. 
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The operator of the axial-vector current is a Dirac pseudovector. For a free particle, it is 
of the nonrelativistic form (// = 0, 1, 2, 3) 



(89) 



where v — p/M is the particle velocity, and &i are Pauli spin matrices. For S — 1, the 
exchange part of the particle-hole interaction is to be taken as 



9 (nn') ^i^ial = -g (nn') , 



(90) 



and 



e+ = v5,o, e = (91) 

Then for a space part of the correlation tensor = 1, 2, 3) we find K'^ — SijKA, where 



and the full vertices are to satisfy the equations 

_ , , If dn' , , 

The temporal component is of the form: 



K (n) r_ (n) - — K (n) T+ (n) 



(92) 



= 1, 



in') r_ (n') - (k in') - 2A (n')) T, (n') 



= 0. 



47r 



qv 



K (n) T° (n) + [k (n) - 2 A (n)) T° (n) 



(93) 
(94) 

(95) 



where the full vertices should be found from the following set of equations 

1 f dn' 



UJ 



n {n!) (n') - (n') (n') 



0, 



UJ 



qv 



-/.(nOT"(n')-(«:(n')-2A(n'))Tf(n') 



= V. 



Mixed space-time components are given by 



^S'(-.q) = P«/| 



UJ 



— Ac(n)T_(n) + («:(n)-2A(n))T+(n) 

V 



dn 

47r 



K (n) r° (n) - (n) (n) 



(96) 
(97) 

(98) 
(99) 
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To obtain a solution in reasonably simple form, we approximate the interaction amplitude 
by its first two harmonics, according g (nn') = gQ + fi'inn'. Then we find the full vertices in 
the form 



l + (7o (2a-r/(l + 5is2r/))/4' 



B2'nscose 

T+ n = 



l-(7or?(l + 52s2ry) /4' 



where 



1 -^70^7(1 + 525^^7) /4' 



B^{uj,q) = ^g,(^l-^g,f3^ , (102) 
52(cu,q) = ^(7i(l-^^?i(/5-27)] • (103) 



Simple algebraic calculations yield the following autocorrelation functions: 

, 2a - ri {1 - Bis^ri) 
1 + go {2a-r]{l + B,s^r]))/4 



and 

T^ij ( \ X 7/(1 + 525^77) 

A y (^^^,0) = -; — „ , . , (105) 

Mixed components and K"^^ are given by the integrals fl95l) and fl99|) . where, according to 
Eqs. ( !23l) . (llOOp and ( llOip . the integrands are odd in cos^. Therefore the mixed polarization 
vanishes: 

K'^icu,q) = K'^{cu,q) = 0. (106) 

Let us consider various limits in the expressions obtained above. For arbitrary tempera- 
ture T > and u, qVp ^ A, according to Eq. (!60!) . we have 

(3 = 2j + s'v-l, (107) 



and 



^^ = J3K A TYTI- (108) 

4 1 _ ^1 (s^r] - 3) /4 
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T 





len the spin-density autocorrelation function fllOSp reproduces the result obtained in Ref. 



where 



A-^(..q) = .,p ^_;W^^^^ . (109) 

and Tj (s) is given by Eq. (!59l) . 

Next we consider the case of time-like momentum transfer when qVp -C A, > 2 A, and 
thus M = s"^ < 1. From Eqs. ([70]), dZI]) and ([7SD, ([75]), we find in this limit 

7(^,r)~^«(cu,r). (Ill) 



For > 0, we obtain 

cuVcu2-4A2 |l + ^o« (a;,T) /2|' 4T 
and 



ImirT(^,q) - r^^. ^^^\^. tanh^, (112) 



where 



lmir^(.,q) ^ -^,%l-S!^!^i^^^l±^lZl^tanh^, (113) 
'^^ *J3^i^ ^3^^2 _ 4^2 |i + (^,a(^,T)/6|' 4r' ^ ^ 

, ^ p & A2 E . A2e(cu-2A) u 

VI. NEUTRINO ENERGY LOSSES CAUSED BY PAIR RECOMBINATION 

As an application of the obtained results we consider the neutrino-pair emission through 
neutral weak currents occurring at the recombination of quasiparticles into the ^5*0 conden- 
sate. The process is kinematically allowed thanks to the existence of a superfluid energy gap 
A, which admits the quasiparticle transitions with time-like momentum transfer q = {u, q), 
as required by the final neutrino pair. 

We consider the total energy which is emitted into neutrino pairs per unit volume and 
time which is given by the following formula (see details, e.g., in Ref. 13l|): 



/G^y^ 2 ImnC-,, (g) Tr (/,/:) d% d^q, 
\2V2J exp(|)-l 2a;i(2vr)3 2.;2(2vr)3' 

where is the Fermi coupling constant, is the neutrino weak current, and ^'^^^y- is 
the retarded weak polarization tensor of the medium. The integration goes over the phase 
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volume of neutrinos and antineutrinos of total energy uj = uji + 002 and total momentum 
q = + q2. The symbol Ylu indicates that summation over the three neutrino types has 
to be performed. 

By inserting J d^q5^^^ {q ^ Qi ^ Q2) = 1 in this equation, and making use of the Lenard's 
integral 

/ ^^^/'^ ~ ~ ^^-^^^ = Y (^'^^^ " ^'^'^^^ ® ® ^""^ ' ^^^^^ 

where (^^jy = diag(l, —1, —1, —1) is the signature tensor, we can write 

e = / du [ dq q2 -ImW^l^^ (g) {q^q, - q^g^J) , (117) 

487?"* Jq Jq exp (yj - 1 

where AC = 3 is the number of neutrino flavors. 

In general, the weak polarization tensor of the medium is a sum of the vector-vector, 

axial-axial, and mixed terms. However, the medium polarization in the vector channel can 

be neglected, because the imaginary part of the longitudinal and transverse polarization 
unctions is proportional to Vp <^ 1, as given by Eqs. fl82l) and fl88l) . (See also Refs. j^, 
^ for details). The mixed axial- vector polarization has to be an antisymmetric tensor, and 

its contraction in Eq. f lll7p with the symmetric tensor q^^q,^ — q^Qf^u vanishes. Thus only 

polarization in the axial channel should be taken into account. 

We then obtain Imll^^gg^^ ~ C\lmK'^, where Ca is the axial weak coupling constant of 

the baryon. Making use of Eqs. (11061) . (11121) . and (I113p . we find 

ImW^ (q) {q,q. - q^g,^) = -^Clp,M*V^^^^Z3§ tanh ^ 

By inserting this into Eq. (11171) and performing integration over dq, we obtain the neutrino 
emissivity in the axial channel, which can be represented in the form 

e = —-GICIKpfM*V^TV / dx- 

15vr5 Jo (e^ + 1) 

.M^ (i + g./i2)^ n I 

' r>'P |1 + Joa (</, z) /2|= 21 |l + g,a (y, z) /6f 



2 
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where y = t^jT and z = y^x^ + y'^. The function a (y, z) is given by 



a {y, z) 



y 



dv y'^ ^ V 

^ TV tanh — 

/y^ - {z^ - v^) 2 

1 z 
= tanh -, 
2 2 



(120) 



Some comments on the approximations done in previous works would be here appropriate. 



In Refs. 



Ill, B, Q, 



17( 1 the calculation of the neutrino emissivity is performed in 



the BCS approximation, i.e. the authors discard Fermi-liquid interactions in a superfluid 
system. The attempt to take into account the particle-hole interactions was undertaken 
recently in Ref. However, though the authors state the important role of the particle- 
hole interactions, their final result for neutrino emissivity contains no Landau parameters 
characterizing this interaction [see Eq.(35) of Ref. [^j. As a matter of fact this means 
that the Fermi-liquid effects have been discarded in this calculation and the result also 
corresponds to the BCS approximation. 

Thus only the BCS limit of our Eq. flllQp can be compared with the previous calculations. 
Setting go = gi = 0, we obtain 



,BCS 



M 



*2 



21) 



+ — ) GlClKvpM^V^T'y^ 



dx- 



\2 ' 



1211 



iQ (e^ + 1) 

where y = /\/T and z = a/x^ + y'^. 

Although this expression reproduces the known BCS result for the neutrino emissivity in 
the axial channel, we recall that the total neutrino emissivity, as given by this formula, is 
suppressed as Vp with respect to the earlier results because the vector channel is practically 
closed. Second term in the brackets was for the first time obtained in Ref. l|. The first term 
is the same as in Ref. ll| . Notice that thj^ term originating from the temporal component 



of the axial-vector current is lost in Ref. 



6|- 



We also do not support the result obtained in Ref. [5|], where one more term is suggested 
due to the mixed space-temporal polarization of the medium. In our calculations, the mixed 
contribution, being odd in cos 6*, vanishes on angle integration: see our Eq. (110611 . This 

n n □ 

agrees with the results obtained in Refs. nj, lUi, 13. 

n n 

The temperature dependence of the energy losses, as obtained in Refs. |5[ , [6| , also is not 
convincing, because the imaginary parts of the polarization functions are calculated for zero 
temperature when no broken Cooper pair exists. The temperature dependence, as given 
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in our Eq 



• (11211) ■ has been repeatedly obtained by many authors before (see, e.g., Refs. 



This dependence follows directly from the kinematics of the reaction and 
statistics of the pair-correlated fermions. 

According to our Eq. flllSp . the imaginary part of the retarded polarization tensor 
substantially depends on the temperature. This dependence may be easily understood in 
the BCS approximation. In this case, 

Imll^'' oc tanh— , 

and (besides the temperature dependence of the energy gap) the temperature-dependent 
factor in the integrand of Eq. flll7p . 

tanh — = ^, (122) 

4T (exp^ + 1)'' ^ ' 

represents the product of occupation numbers in the initial state of two recombining quasi- 

particles. Indeed, the dominant contribution to the phase integral enters from the quasipar- 

ticle momenta near the Fermi surface. As the neutrino-pair momentum q ~ Tc -C pi?, one 

can neglect q in the momentum conservation 6 function, thus obtaining p' = —p. After this 

simplification, the energies of initial quasiparticles are Ep' — Ep — /2 . 

VII. NUMERICAL EVALUATION 

In Eq. flligp . the temperature dependence of the emissivity enters by means of parameter 

A(r) A(0)A(r) 

with T = T/Tc, where T^ is the superfiuid transition temperature. For a singlet-state pairing 



A (0) /Tc = 1. 76 (see, e.g., Ref. [l4|), therefore the function y depends on the dimensionless 
temperature r only. Thus, the emissivity in Eq. flll9p . in the standard physical units, can 
be written as 
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where Mp is the bare proton mass, C\ = ~ 1.6 (for neutrons) , and the function F (r) is 
defined as 

F{T) = T'y / dx -2 

/gV^ + 1 j 

^ / M*^ (1+W12)' 11 1 ^ .^25) 

\M^\l + g,a{x,y)/2\' 21 \l + g,a{x,y) ) 

The function a (x, y) can be recast as 

+ .^^I^tanhV^?±2. (126) 
4 xa/x^ + 2 

In numerical estimates, we use the fit expression of the energy gap dependence on the 
temperature (see, e.g., Ref. |ll|): 



Unfortunately, the Landau parameters g^, gi are poorly known up to now. These are 
known to depend on the baryon density and could be of the order of unity |l^, [l9]. Ex- 
tracted from nuclear data, g^ = 1.5, while gi is unknown ^^]. In our estimate, we use three 
different combinations of these parameters. The result of numerical evaluation is shown in 
FIG. [H where we compare the energy losses according Eq. 0124p with the BCS expression 
(11211) . which can be cast in the same form as Eq. fll24p but with the function F (r) replaced 
by 

Fbcs (r) = rV H dx (^^')' (128) 
^0 |gV^^+f^ + 1 ) 

This function is represented by the lowest curve. The upper curves represent the ratio 
F (r) / -Fbcs (''') for three different combinations of the Landau parameters. 



VIII. SUMMARY AND CONCLUSION 



In this paper, we have investigated the Fermi-liquid effects in the neutrino emission 
at the pair recombination of thermal excitations in a superfluid crust of neutron stars. 
For this purpose, we have calculated the weak response functions of superfluid fermion 
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FIG. 1: The temperature dependence of neutrino energy losses. Lowest curve ~ the function FbcSj 
as given by Eq. (jl28p . Upper curves - the ratio -F/-Fbcs for three different combinations of Landau 
parameters go,gi shown near the curves. 

system at finite temperatures while taking into account the particle-hole interactions near 
the Fermi surface. For the calculation, we used Leggett approach to strongly interacting 
Fermi liquid with pairing. In the case qVp -C A, typical for the weak processes in the 
nonrelativistic baryon matter of neutron stars, we have derived the response functions valid 
at finite temperature and for arbitrary transferred energy uj ^ A. Our general expressions, 
as given by Eqs. (!55|) . (!88|) . fll04p . and (llOSp . naturally reproduce the well-known results 
jol, [lO| obtained for the case of small transferred energy, ^ A, as well as the response 
functions obtained for arbitrary uj in the BCS approximation 4|. 

In the kinematic domain uj > 2 A and q < cu, we have carefully calculated the imaginary 
part of the response functions up to the necessary accuracy, what allows us to evaluate 
the neutrino energy losses caused by the pair recombination while taking into account the 
Fermi-liquid effects. 
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In the vector channel, we found that the spherical harmonic of the particle-hole interac- 
tion does not affect the imaginary parts of polarization functions in the time-like domain. 
The imaginary part of both the longitudinal and transverse polarization functions are pro- 
portional to Vp, and thus the particle-hole interactions are not able to increase substantially 
the intensity of neutrino-pair emission through the vector channel. 

The imaginary part of the axial polarization is suppressed as Vp, therefore the dominating 



neutrino emission occurs in the axial channel. We do not support the statement of Ref. 



5| 



that the particle-hole interactions can be ignored [see the discussion after Eq. (33) of Ref. 
jsl). Our analytic expression (11241) and numerical estimates demonstrate the important role 



of the Fermi liquid effects in the considered process. 



IX. HERE 

Discarding the particle- hole interactions means that the result obtained in Ref. [sl, as 
a matter of fact, corresponds to the BCS approximation. This approximation has been 
used before by several authors. Therefore for comparison, we consider the BCS limit of our 
Eq. (11191) which can be obtained by putting go = di = 0. The detailed analysis of some 
controversial results of different authors can be found at the end of Sec. VI. 

For completeness, it is helpful to discuss additionally the case, when the quasiparticles 
carry an electric charge. Though the direct neutrino interaction with recombining protons 
is screened by the proton background ^, the proton quantum transitions can excite back- 
ground electrons, thus inducing the neutrino-pair emission by the electron plasma. This 
effect has been already studied in Refs. 12|, 13|; therefore, we only briefly revisit this 
problem in the light of modern theory to understand whether the plasma effects can lead 
to noticeable neutrino energy losses through the vector channel. For the sake of simplicity 
we consider a degenerate plasma consisting of nonrelativistic superfiuid protons and rela- 
tivistic electrons. As found in Refs. [l^, [3], the role of the electron background, in this 
case, consists of the effective renormalization of the proton vector weak coupling constant, 
c'y jl c~y ■ Thus we find that the electron background strongly increases the effective 
proton vector weak coupling with the neutrino field, {^c"y l(^y^ ^ 576. However, this huge 
factor should not mislead the reader, because it arises only as a result of a very small proton 
coupling constant, c'y <^ c~y . Since the degenerate electron plasma can be considered in the 
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collisionless approximation, the imaginary part of the medium polarization arises from the 
proton pair recombination and therefore is proportional to Vp, where Vp -C 1 is the Fermi 
velocity of protons. Thus the neutrino emission through the vector channel is suppressed by 
a small factor Vp and may be ignored in comparison with the dominating neutrino radiation 
in the axial channel, where the neutrino energy losses are suppressed as Vp. 

We now return to the Fermi-liquid effects incorporated in Eq. flllQp . The magnitudes 
of the Landau paramete rs Q o, Qi are poorly known and depend on the baryon density. By 
modern estimates [l8|, [l9|, these could be of the order of unity. Thus the Fermi-liquid 



effects can notably modify the emissivity dependence on the temperature and the matter 
density as compared to that found in the BCS approximation. This, however, cannot change 
the main conclusion that the dominating contribution to the neutron and proton emissivity 
comes from the axial channel of weak interactions . This means that the neutrino energy 
losses are to be suppressed as compared to that of Ref. [ll by a factor of Vp. This could serve 
by a natural explanation of the observed superburst ignition discussed in the Introduction. 
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